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HEAT SURGERY SYSTEM MONITORED BY REAL- 
TIME MAGNETIC RESONANCE PROFILING 

CROSS REFERENCE TO RELA TED APPLICATIONS 

This application is related to U.S. Patent Applications 
5 "Magnetic Resonance Guided Pocussed Ultrasound Surgery" by 
Harvey Cline, Robert Ettinger, Kenneth Rohling and Ronald 
Watkins, Ser. No. 07/854,040 filed 3/19/92, "Magnetic Resonance 
Surgery Using Heat Waves Produced with A Laser Fiber or 
Focussed Ultrasound" by Harvey E. Cline and Thomas R. Anthony 

10 et aL Ser. No. 07/751,259 filed 8/29/91 and U.S. Patent 5,133,357 
issued July 28, 1992 "Quantitative Measurement of Blood Flow 
Using Cylindrically Localized Fourier Velocity Encoding" by 
Charles L Dumoulin, Christopher J. Hardy, Steven P. Souza, and 
Steven A. Ash all assigned to the present assignee and hereby 

15 incorporated by reference. 

BACKGROUND OF TffE INVENTION 

The present invention relates to surgery performed by local 
heating guided by magnetic resonance (MR) imaging, and more 
particularly to a method and apparatus for performing surgery by 
20 pulsed local heating guided by real-time temperature-sensitive MR 
imaging. 

Conventional MR imaging provides an operator, such as a 
surgeon or a radiologist, with internal views of a patient's 
anatomy. MR imaging provides excellent contrast between 
25 different tissues and is useful in planning surgical procedures. 
Many tissues in a patient such as a tumor, are much more visible 
in an MR image than as seen in actual surgery. The tumor can 
also be obscured by blood during surgery further reducing 
visibility. 

30 Experiments on living tissue show that a heated zone above a 
critical temperature destroys tissue. This zone increases in size 
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with time as the heat is applied to reach a steady state of both 
temperature and heat flow. If the maximum temperature is 
limited, then an area heated to a temperature exceeding the 
critical temperature causing destruction of tissue, approaches a 
5 predetermined size. 

Researchers at Brigham and Womens Hospital, Boston, MA 
have proposed treatment of deep lying tumors by laser surgery, as 
described in "MR Imaging of Laser-Tissue Interactions 1 *, by FA. 
Jolesz, A.R. Bleire, P. Jakob, P.W. Ruenzel, K. Hutti, G.J. Jako, 
10 Radiology 168:249 (1989). 

In laser surgery of the brain, a small burr hole is drilled in the 
skull and a hollow needle containing an optical fiber is then 
inserted into the tumor. Optical energy is passed through the 
optical fiber and into the tumor heating the tumor. Lasers guided 
15 through fiber optics are potentially very useful for surgery in the 
brain, since they allow access (through biopsy needles) to deeply 
buried pathology with minimal disruption of intervening tissue. 
The laser destroys the pathological tissue through localized 
heating. 

20 A view of the heated region is provided with the use of MR 
temperature-sensitive pulse sequences. One known MR 
temperature-sensitive pulse sequence is described in U.S. Patent 
4,914,608 "In-vivo Method for Dete rmining and T m a gin g 
Temperature of an Object/Patient from Diffusion Coefficients 

25 Ob taine d by Nuclear Magnetic Resonance*\by Denis LeBihan, Jose 
Delannoy, and Ronald L. Levin issued April 3, 1990. This sequence 
causes full MR images to be created. It is, therefore, relatively 
slow, time-consuming and not capable of monitoring quickly 
changing temperatures. 

30 In focussed ultrasound surgery, acoustic energy is 

concentrated at a focal point within the tumor. The focussed 
acoustic energy causes heating of tissue at the focal point 
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Tumors have been selectively destroyed in cancer patients 
using focussed ultrasound to heat tissue within the patient at the 
University of Arizona as reported in "Effects of Physical 
Parameters on High Temperature Ultrasound Hyperthermia" by 
5 B. E. Billard, K. Hynynen and Robert. B. Roemer, Ultrasound in 
Med. & Biol. Vol. 16, No. 4, pp. 409-420, 1990. Billard et al. disclose 
that the control of heat is improved by using short heating pulses 
where the effect of blood perfusion is negligible. However, since 
they do not image the temperature distribution, it is difficult to hit 
10 small, deep laying targets. 

It would be beneficial to accurately image, in real time, 
localized heating and selectively destroy a specific tissue with 
minimal invasiveness without harming surrounding healthy 
tissue. 

15 OBJECTS OF THE INVENTION 

It is an object of the present invention to selectively heat tissue 
of a patient with a small degree of invasiveness and accurately 
monitor the heated tissue in real time with magnetic resonance 
imaging. 

20 It is another object of the present invention to accurately destroy 
selected tissue without harming surrounding tissue with a non- 
invasive procedure employing real-time magnetic resonance 
temperature-sensitive imaging. 

SUMMARY OF THE INVENTION 

25 In accordance with the preferred embodiment, a pulsed heat 
source concentrates heat at an application point. A positioning 
means, operated by an operator, such as a radiologist or surgeon, 
positions the application point to selectively heat specific tissue of a 
patient with minimal invasiveness. A magnetic resonance (MR) 

30 imaging device employs a real-time temperature-sensitive MR 
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pulse sequence to create real-time temperature-sensitive images of 
the heated tissue which are provided to the operator performing 
the procedure. The operator then interactively adjusts the amount 
of heat concentrated at the application point and the position of the 
5 application point to selectively destroy the specific tissue. 

The temperature-sensitive real-time MB pulse sequence causes 
an elongated region of magnetization to be excited within a specific 
tissue of the patient, employing a cylindrical excitation, or 
excitation of the intersection of two imaginary slices. The excited 

10 tissue emits an MR response signal which is temperature encoded 
by with a Ti-weighted readout gradient, or diffusion of magnetic 
field gradients. A one-dimensional Fourier transformation (1DFT) 
is applied to the MB response signal to generate a magnetization 
profile. This profile indicates the temperature of tissue along the 

15 length of the excitation region. The excitation region is oriented to 
intersect tissue at the application point, thereby permitting real- 
time monitoring of heat destruction of the tissue. 

BRTEF DESCR TPTTON OF THE DRAWINGS 

The features of the invention believed to be novel are set forth 
20 with particularity in the appended claims. The invention itself, 
however, both as to organization and method of operation, together 
with further objects, may best be understood by reference to the 
following description taken in conjunction with the accompanying 
drawing in which: 

25 Pig. 1 is a block diagram of a magnetic resonance (MR) pulsed 
heat system according to the present invention. 

Fig. 2 is a partial illustration of a first embodiment of the present 
invention employing focussed ultrasound as a pulsed heat 
means with a patient positioned for surgery within the bore 
30 of the magnets . 



WO 94/23308 



PCT/US94/02367 



-5- 

Fig. 3 is a partial illustration of a second embodiment of the 

present invention employing a laser source and fiber optics 
as a pulsed heat means with a patient positioned for surgery 
within the bore of the magnets. 

5 Fig. 4 is a representation of a bulk magnetization vector M, here 
equal to an equilibrium magnetization vector Mo arising 
from unpaired resonant nuclei ("spins") within the tissue of 
a patient in a static magnetic field, shown in a reference 
frame XTZ* rotating about the static magnetic field axis at 
10 a frequency equal to the precession frequency of the spins. 

Fig. 5 is a representation of the bulk magnetization M of Fig. 4, 

also in reference frame XTTZ\ after an MR excitation pulse 
is applied to tip it into the transverse plane XT, so that now 
the transverse component of magnetization MxY = M. 

15 Fig. 6 is a graphical representation of a resultant magnetic-field 
gradient vector 6 r> and its trajectory over time during 
cylindrical excitation. 

Fig. 7 is a cylindrical-selection TV weighted real-time temperature- 
sensitive MR pulse sequence according to the present 
20 invention. 

Fig. 8a is a graphical illustration of the shape of a region of tissue 
excited by cylindrical-selection pulses. 

Fig. 8b is a illustration of a cross-section of transverse 
magnetization of the region of Fig. 8a. 

25 Fig. 9 is an illustration of an X component of the bulk 

magnetization vector M in the rotating frame of Fig. 5 
during an applied magnetic field gradient. 



Fig. 10 is a timing diagram of an intersecting slice selection and 
diffusion gradient encoded real-time temperature-sensitive 
MR pulse sequence according to the present invention. 

Fig. 11 is an illustration of the intersection of two imaginary slices 
according to the present invention. 

Fig. 12 is a cylindrical-selection diffusion gradient encoded real- 
time temperature-sensitive MR pulse sequence according to 
the present invention. 

Fig. 13a is an axial MR image of an agarose gel phantom 
surrounding a tube of cold water. 

Fig. 13b is an axial MR image of the same agarose gel phantom of 
Fig. 13a surrounding a tube of hot water. 

Fig. 13c and 13d is an image of a series of temperature vs. position 
graphs, temperature-sensitive profiles* of the phantoms of 
Figs. 13a, 13b, respectively. 

Fig. 13e is a temperature vs. position graph of a single 
temperature-sensitive profile of Fig. 13c. 

Fig. 13f is a temperature vs. position graph of a single 
temperature-sensitive profile of Fig. 13d. 

Fig. 14a is a position vs. time temperature-sensitive image 

acquired with the present invention of the focal point of a 
focussed ultrasound pulsed heat means of the present 
invention. 



Fig. 14b is a position vs. time temperature-sensitive image being a 
continuation of Fig. 14a. 
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DETAILED DESCRIPTION OF THE INVENTION 

A specific tissue, such as tumor tissue, in a patient may be 
selectively destroyed by localized heating without affecting the 
surrounding healthy tissue. A method of controlling the size of a 
5 region being heated by pulsing a heat source is disclosed in 

"Thermal Diflusity of Isotopically Enriched 12 C Diamond" by T. R. 
Anthony, W. P. Banholzer, and J. F. Fleischer Phys. Rev. B Vol. 
42, No. 2 July 15, 1990, hereby incorporated by reference. 
Similarly, in the present invention heat is applied to create a 

10 heated region in the tumor tissue in a pulsed or oscillating 
fashion. This oscillation creates a heat wave at an application 
point. The pulsed heat may be produced by a focus sed ultrasound 
transducer or a laser and fiber optics. (A more detailed 
description of the heating profile of the ultrasound transducer is 

15 included in the aforementioned patent application Ser. No. 
07/854,040). The heated region is monitored with a magnetic 
resonance (MR) imaging system employing a temperature* 
sensitive real-time MR pulse sequence. Images of internal 
structures of the patient and the heated regions are supplied to an 

20 operator allowing the operator to interactively adjust the position 
and size of the heated region. 

A block diagram of the MR pulsed heat system of the present 
invention is shown in Fig. 1. An MR imaging means 10 comprises 
a static field magnet 260, 270, a gradient means 40, an RF 

25 transmitter 50, a receiver means 90 and a computation means 95 
which rapidly acquire images of a patient 200. Static field magnet 
260, 270 provides a static magnetic field Bo over patient 200 required 
for MR imaging. Gradient amplifier 40 and RF transmitter 50 
supply the necessary magnetic field gradients and RF radiation, 

30 respectively, according to MR pulse sequences over patient 200 
causing portions of patient to emit an MR response signal. Raw 
MR response signals are sensed by receiver 90 and passed to a 
computation means 95 which computes an MR image. 



WO 94/23308 



PCT/US94/02367 



Computation means 95 may then display images to an operator on 
a display means 110. A control means 100 receives MR images 
from computation means 95. The operator interfaces with control 
means 100 to compute a path from pulsed heat means to a desired 
5 location within patient 200 which avoids obstructions such as bone 
and air spaces within patient 200. In the case of brain tumors, a 
conventional MR image is acquired to locate the tumor and plan a 
safe trajectory between an entry point and target point within the 
tumor. 

10 Control means 100 then actuates positioning means 70 to 
position pulsed heat means 80. MR imaging means 10 employs 
real-time temperature-sensitive pulse sequences to rapidly acquire 
temperature-sensitive images of patient 200 which may optionally 
be superimposed on another medical diagnostic image (which may 

15 be the conventional MR image). Since both the internal structures 
and heated regions are imaged, the operator can accurately 
position the heated region to correspond to a desired internal 
structure, and accurately heat specific tissue. 

A first embodiment of the MR pulsed heat system according to 
20 the present invention employing a focus sed ultrasound transducer 
as the pulsed heat means 80 is shown in Fig 2. Patient 200 is 
placed on a table 310 designed to accommodate a focus sed 
ultrasound transducer 330 in an ultrasound conducting liquid 
bath 320. The ultrasound conducting liquid chosen is one that will 
25 conduct ultrasonic energy with little attenuation. The ultrasound 
transducer 330 is moved inside the bore of static field magnets 260, 
270 by positioning means 70 to focus on different locations within 
patient 200. A path is computed by control means 100 under the 
direction of the operator from a set of images of the patient 
30 avoi ding bone or air in the path of the ultrasound beam. The focal 
point of ultrasound transducer 330 is positioned by positioning 
means 70 onto a tissue desired to be heated or destroyed, such as 
tumor 280 and pulsed to selectively heat tumor 280. The 
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ultrasound transducer is interactively positioned by the operator 
while the operator views cross sectional temperature-sensitive 
images. (A more detailed description of the positioner is included 
in the aforementioned patent application Ser. No. 07/854,040). 

5 - Another embodiment employing a laser and optical fiber as 
pulsed heat means 80 is shown in Fig. 3. Patient 200 lies on a table 
210 that moves into the bore of the two part static field magnet 260, 
270. A laser fiber 230 is inserted into the patient with a hollow 
needle 240 which may optionally be guided by a mechanical 

10 positioning device 250 driven by the operator through control 
means 100 of Pig. 1. The trajectory of hollow needle 240 may be 
computed from conventional MR images of the patient which does 
not intersect critical anatomy such as large blood vessels. Heat is 
applied to tumor tissue 280 by periodically pulsing the laser 

15 through laser fiber 230 (i.e., a fiber optic material) to selectively 
destroy tumor 280 while the operator views a real-time 
temperature-sensitive MR imag e. More than one needle may be 
required to remove an irregularly shaped tumor. 

Free impaired spinning protons ("spins") in the nucleus of a 
20 molecule of a specific tissue, normally hydrogen nuclei, align 
themselves in a magnetic field Bo such that their axes precess 
about the magnetic field* Each spin is a small magnetic dipole. 
The net sum of the population of dipoles results in a bulk 
ma gne tization M which is aligned with static magnetic field Bo, 
25 shown in Fig. 4 in a reference frame XTT2T rotating about the 
static magnetic field axis at a frequency equal to the precession of 
the spins . A radio frequency (RF) pulse applied in the presence of 
a magnetic field of predefined strength causes excitation or 
resonance of the spins, causing M to tip into the transverse plane 
30 XT, thereby increasing transverse magnetization vector M x y of 
the spins, as shown Fig. 5. 
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It is possible, by choosing the strengths of the RF pulse and 
magnetic field gradients, to selectively choose spins for excitation. 
Spatial localization of spins may be performed by choosing an RF 
pulse and magnetic field gradients to select spins of a specific 
5 desired area of the patient to be imaged. 

The present invention achieves spatial excitation by either 
cylindrical-selection, or intersecting slice selection. In a first 
embodiment, a cylindrical region of the patient is excited or 
selected with RF pulse 31 and gradient pulses 33, 35 of "X" and "Y" 

10 directions, respectively, as shown in Fig. 7. The "XT and T 

directions pertain to two orthogonal directions, both perpendicular 
to the axis ("Z M direction) of the cylindrical excited region referred 
to as the logical coordinates. Logical coordinates may be converted 
to physical coordinates which are stationary with respect to 

15 physiral objects, such as the MR imaging apparatus, by a 

conventional coordinate transformation. RF pulse 31 changes 
amplitude over time as the G r and G y orthogonal magnetic field 
gradients 33 and 35 create a resultant gradient field vector Gr as 
shown in Fig. 6. The tip of vector Gr traces a spiral path. The G x 

20 and Gy orthogonal gradient waveforms 33, 35 and RF pulse 31 
applied simultaneously, result in excitation of a cylindrical 
element of the patient. RF pulse 31 of Fig. 7 is chosen to be the 
weighted two-dimensional Fourier transform of the desired 
excitation profile. The method of determining gradient waveforms 

25 and RF pulses to excite desired shapes is disclosed in U.S. patent 
application "Sampling-Ring Saturation Pulse for Two- 
Dimensional Magnetic Resonance Selective Excitation" by 
Christopher J. Hardy, and Edward J. Nieters Ser. No. 07/928,390 
filed 8/12/92, assigned to the present assignee and incorporated by 

30 reference. 



After application of RF pulse 31 and magnetic field gradients 
33, 35, excitation cylinder 81 has a transverse magnetization as 
indicated in Figs. 8a and 8b. The center of cylinder 81 is at 
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displacement x* At displacement xi of Fig. 8a, the transverse 
magnetization M x y increases to m f as shown in Fig. 8b. 
Transverse magnetization M x y remains substantially constant 
at m across the diameter of cylinder 81, which drops off at radius 
5 X2 to zero. 

After the cylinder 81 is excited, transverse magnetization M x y 
of Fig. 5 processes in the transverse plane. As shown in Fig. 7, a 
first portion 37a of readout gradient pulse 37 applied at time ti 
causes spins to become partially dephased such that spins are 

10 pointing in different directions within the transverse plane at any 
instant in time. A second portion 37b, essentially having the same 
size shaded area (the pulse duration multiplied by the gradient 
amplitude) of readout gradient pulse 37 applied at time t2» causes 
spins to become rephased and peak at time ts creating a half-echo 

15 MR response signal 39. A last portion 37c of magnetic gradient 
pulse 37 causes the MR response signal to be encoded along the 
direction of gradient 37c A longitudinal component Mz of 
magnetization, aligned along the static magnetic field Bo, is 
related to a spin-lattice relaxation time Ti of the tissue by the 

20 following equation: 



M,= 



' Mq for the first MR response signal 

Mq^I- 6 T> j for subsequent MR response signals 



[1] 



where Mo is the equilibrium longitudinal magnetization of the 
spins, and Tr is the time before the next repetition of the pulse 
sequence. The longitudinal magnetization Mz grows toward the 

2 5 equilibrium value Mo with a time constant Ti, which is affected by 
temperature. Ti is increased in tissue by approximately 1% for 
each 1° C change in temperature, due to thermal influences on 
spin-lattice interactions. RF excitation pulse 31 v of subsequent 
repetitions tips Mz into the transverse plane, producing an MR 

30 response signal M s i g which is indicative f temperature of the 
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tissue in the excited r gion, A one-dimensional Fouri r 
Transform (1DFT) is applied to the MR response signal to 
determine the amplitude of the MR response signal for each 
frequency, which is mapped to each location along the readout 
5 gradient, thereby producing a temperature-sensitive profile along 
the direction of the readout gradient 

This entire pulse sequence is repeated after time Tr, starting 
with RF pulse 31' and gradient pulses 33' and 35'. 

In a second embodiment of the invention, diffusion sensitivity is 
10 encoded. Fig. 9 is an amplitude vs. time diagram of the "X m 

component My of the bulk magnetization M at a single location, 
being a component of the total transverse magnetization vector 
M x y of Fig. 5 within a tissue. The vector M rotates at a specific 
constant frequency during period 54 related to the magnetic field 
15 experienced by the spins. A magnetic field gradient 57 having an 
amplitude G is applied during period 55, changing the magnetic 
field experienced by the spins, thereby changing the frequency of 
the rotation of the bulk magnetization vector M. During period 56 
of Fig. 9, M reverts back to its original frequency, but the phase has 
20 been increased. This is known as phase evolution. 

The relative phase evolution of the transverse magnetization 
vector M x y of Fig. 5 at each point in space is determined by the 
deviation in magnetic field from value BO at that point (which is in 
turn determined by the amplitude of the applied gradient) and the 
25 time during which it is applied. The relative phase shift of the 

transverse magnetization arising from stationary spins is directly 
proportional to the lobe area under magnetic field gradient pulse 
57. 

By applying a gradient to a group of moving spins, causing 
30 phase evolution in one sense, followed by application of another 
gradient, causing phase evolution in an opposite sense, a net 
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phase evolution may be encoded indicating the motion of the spins. 
This phase shift, ^(motion), may be described as: 

4><motion)=YVTAg [2] 

where y is the gyromagnetic ratio specific to a given element, V is a 
5 component of the spins* velocity parallel to the direction of the 
applied gradients, T is the time between the centers of the 
magnetic gradient pulses applied along the line in which flow is to 
be measured, and A g is the area (gradient strength multiplied by 
duration of its application). Equation [2] ignores phase shifts 
10 which can arise from higher orders of motion such as 
acceleration, and change in acceleration, etc. 

Since diffusion is motion of spins, and temperature has a direct 
relation to the rate of diffusion, diffusion encoding of the type 
described above may be employed in measuring temperature. 

15 Fig. 10 illustrates a second embodiment of the real-time 

temperature-sensitive MR pulse sequence of the present invention 
employing intersecting slice selection and diffusion gradient 
encoding. An RF sine(xyx ("sine") pulse 41 is applied 
simultaneously with an X magnetic field gradient pulse 42 

20 causing spins to nutate until they are oriented in the transverse 
plane (a 90° pulse). Cross-hatched portions of pulse 42 are of equal 
areas above and below the axis. Sine pulse 41 transforms in a 
frequency domain into a narrow band square region from 
encompassing a small frequency range. Application of RF sine 

25 pulse 41 with gradient pulse 42 has the effect of exciting a narrow 
slice of spins of the patient in the logical Y^Z plane processing at 
corresponding frequencies. Diffusion gradients 43x, 43y, 43z are 
applied to cause phase evolution in a forward sense. 

A RF refocussing pulse 44 acts to reverse the relative phase of 
30 each of the spins. If a spin has gained 80° of forward phase 

evolution during application of the diffusion gradient, then after 
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the application of RF refocussing pulse 44, the same spin now has 
a -80° phase, or $ = - 

A Y magnetic field gradient is applied simultaneously with RF 
refocussing pulse 44 to select a region, that region being the 
5 intersection of a slice in the X,Z plane and the first selected slice of 
tissue in the Y,Z plane shown as rod-shaped excitation region 82 in 
Fig. 11. 

A second set of diffusion gradients 46x, 46y, 46z is applied along 
the X, Y, Z directions respectively causing another forward phase 

10 evolution. Since RF refocussing pulse 44 has flipped each phase to 
its negative value, a forward phase evolution will cause stationary 
spins to become rephased at time t*. Moving spins will be phase 
encoded by the amount of diffusion in the direction of the net 
gradient, indicative of the temperature. A readout magnetic field 

15 gradient pulse 47 oriented along the length of rod 82 of Fig. 11 (Z 
axis) encodes a full echo MR response signal 58 radiated by rod 82 
which peaks at time t$. Cross-hatched portions of pulse 47 have 
equal areas under the curve. This results in a temperature* 
sensitive profile through rod 82 which may be interactively 

20 reoriented. 

Fig. 12 illustrates a third embodiment of the real-time 
temperature-sensitive MR pulse sequence employing cylindrical 
excitation and diffusion gradient temperature coding. Two 
orthogonal time varying magnetic field gradients 53, 55 and an RF 
25 excitation pulse 51 are applied simultaneously to excite a cylinder 
of tissue within the patient similar to the excitation of tissue 
described in conjunction with Fig. 7. 

Simultaneous diffusion gradients are applied in the X,Y,Z, 
directions 53x, 53y, 53z respectively. An RF refocussing pulse 52 is 
30 applied. A second set of diffusion gradients 56x, 56y, 56z are then 
applied along the X,Y,Z, directions, respectively. An MR response 
signal is then acquired in the presence of a readout gradient 57 
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oriented along the length of the cylinder. A one-dimensional 
Fourier Transform (1DFD is applied to the MR response signal 
resulting in a temperature-sensitive profile along the direction of 
readout magnetic field gradient from magnetic field gradient 
5 pulse 57. Cross-hatched areas of pulse 57 have equal area above 
and below the axis. In this manner a temperature-sensitive 
profile through the specific tissue along the axis of the cylindrical 
excitation region may be recorded once every 50 msec or less. 

The cylinder axis can be oriented in an arbitrary direction by 
10 proper mixing of the X,Y, and Z magnetic field gradients. In 

addition, the cylinder can be offset from the center of the magnet to 
any specified location by frequency modulation of the RF waveform 
of the excitation pulse as described in "Off-Axis Spatial 
Localization with Frequency Modulated Nuclear Magnetic 
15 Resonance Rotating (p) Pulses" by C. J. Hardy, P. A. Bottomley, P. 
B. Roemer, J. Appl. Phys. 63, 4741 (1988), which is hereby 
incorporated by reference. 

Experimental Results 

Experiments indicate that when a laser which emits in the 
20 infrared range of approximately 10W power is used to irradiate 
brain tissue for several minutes, a number of distinct regions of 
tissue damage arise. Cells nearest the end of the fiber optic, which 
experience temperatures >160°F, suffer instant damage and death. 
Those in the immediately surrounding regions, which reach 
25 temperatures > 140°, are also damaged irreparably, but with cell 
death occurring more slowly. In both of these cases the changes 
are readily seen as regions of decreased intensity in an MR image. 
Cells that are heated to temperatures somewhat lower than 140° 
also experience chang es leading to reduced intensity in the MR 
30 image, but in this case the changes are reversible, with recovery to 
a normal MR intensity as the cell cools back down. It is obviously 
very important that the operator is able to distinguish between the 
second and third types of change. Nearby vasculature in the brain 
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acts as a sort of radiator, providing differing degrees of cooling to 
different regions depending on the specific geometry. Therefore, it 
is important to have a means of real-time monitoring of the 
healing process. 

5 A surgical procedure would involve irradiating the specific 
tissue (e.g. tumor) until the region of reduced intensity in the 
temperature-sensitive profile roughly matched the boundaries of 
the tumor. The cylindrical or rod-shaped excitation region could 
be reoriented interactively to probe the tumor in all three 

10 dimensions. A delay allowing recovery to thermal equilibrium, 
again monitored by the MR temperature-sensitive real-time MR 
pulse sequence, would reveal those areas on the border which were 
not irreversibly damaged. The tumor would then undergo 
additional irradiation, perhaps with several iterations of 

15 irradiation and recovery, until the border regions were destroyed 
as well. In this way the amount of damage to surrounding normal 
tissue would be minimized. 

The present invention was employed on a phantom for real- 
time monitoring of laser surgery on a standard SIGNA® MR 

20 diagnostic system* The pulse sequence employed is that of Fig. 7. 
Figures 13a and 13b are axial MR images of the same agarose 
phantom surrounding a tube of water but acquired at two different 
temperatures. The water is circulated from a heat bath using a 
peristaltic pump, and maintained at constant temperature. The 

25 two dark spots in Fig. 13a show the locations of two thermocouple 
wires used to measure the temperature in the agarose. 
Temperature profiling of the present invention was performed in a 
region just superior to the region containing the wires. 

Figures 13c and 13d are temperature-sensitive profiles taken 
30 through the center of the phantom shown in Figs. 13a, 13b, 

respectively, as a function of time, with time running from top to 
bottom. Data for each temperature-sensitiv profile was acquired 
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in 100 msec, with approximately 256 temperature-sensitive profiles • 
taken. Pig. 13c was acquired at room temperature (22.3 C). In Pig. 
13d, the thermocouple just outside the inner tube read 41.6 C and 
the outer thermocouple read 25.9 C. 

5 Figures 13e and 13f are individual temperature vs. position 
profiles taken from Figs. 13c and 13d, respectively. The most 
striking effect of elevated temperature is the drastic drop in signal 
in the inner tube, presumably caused by convection in the water. 
In addition, there is signal loss in the agarose which varies with 
10 distance from the hot inner tube. 

Figures 14a and 14b illustrate the focussed ultrasound 
embodiment of the present invention used to heat and monitor a 
gel phantom. Here, the excitation cylinder has been oriented to 
pass through the ultrasonic transducer's focal point, orthogonal to 

15 the ultrasound beam. The time axis runs from left to right and the 
vertical axis represents the position along the excitation cylinder. 
The ultrasound beam was turned on and off as indicated in Figure 
14a. Rapid heating is evident at the focal point, followed by 
diffusion of the heat and then further diffusion and cooling after 

20 the beam is turned off as shown in Fig. 14b. Fig. 14b was acquired 

immediately after Fig. 14a. 

i 

By employing MR temperature-sensitive pulse sequences 
according to the present invention, temperature can be monitored 
with time resolution as low as 30 msec. The present invention, 
25 therefore, is a viable means of performing heat surgery accurately 
and effectively. 

While several presently preferred embodiments of the invention 
have been described in detail herein, many modifications and 
variations will now become apparent to those skilled in the art. It 
30 is, therefore, to be understood that the appended claims are 
intended to cover all such modifications and variations as fall 
within the true spirit of the invention. 
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WHAT IS CLAIMED IS: 

1. A magnetic resonance (MR) pulsed heat system for allowing an 
operator to selectively heat tissue within a patient comprising: 

a) pulsed heat-producing means adapted for concentrating 
energy at an application point; 
5 b) positioning means for positioning the application point of the 
pulsed heat-producing means in a specific tissue within the 
patient so as to create a heated region within the specified 
tissue; 

c) operator responsive control means for enabling said operator to 
1 0 control the positioning means; 

d) MR imaging means for providing real-time temperature- 
sensitive MR images of the heated region; and 

e) display means for displaying the temperature-sensitive images 
to said operator. 

2. The MR pulsed heat system of claim 1 wherein the MR imaging 
means comprises: 

a) radiofrequency (RF) transmitter for transmitting with a 
repetition time Tr a plurality of RF excitation pulses into said 

5 patient; 

b) a gradient means for: 

1. producing a plurality of time-varying magnetic field 
gradients in an 'IT direction and a plurality of time- varying 
magnetic field gradients inaT direction orthogonal to the 

10 T direction, both gradients applied simultaneously with 

each RF excitation pulse thereby exciting transverse * 
magnetization of an elongated excitation region in said 
patient, 

2. applying with repetition time Tr a plurality of dephasing 

15 magnetic field gradients of a predetermined duration having 

a predetermined amplitude of a first polarity, and 
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3. applying with repetition time Tr a plurality of readout 
magnetic field gradients each having a second polarity 
opposite the first polarity, and having a product of duration 
20 multiplied by amplitude substantially larger than a product 

of the amplitude and duration of the dephasing magnetic 
field gradient; 

c) receiver means for receiving with repetition time Tr a plurality 
of MR response signals from the elongated excitation region 

23 simultaneously with the application of each readout gradient; 
and 

d) computation means for computing a temperature-sensitive 
profile from each MR response signal from the receiver 
means. 

3. The MR pulsed heat system of claim 1 wherein the MR inn* ping 
means comprises: 

a) a radiofrequency (RF) transmitter for transmitting an RF 
sineCxVx "sine" selection pulse and a 180° RF sine refocussing 

5 pulse into said patient; 

b) a gradient means for. 

1. creating a slice selection magnetic field gradient along an 
T direction simultaneously with application of the RF sine 
selection pulse to cause excitation of a first slice of tissue 

10 within said patient, 

2. applying a first set of diffusion gradients in the "X" direction, 
and in the "V and "Z? directions orthogonal to the T 
direction after selection of the first slice and before 
application of the 180° RF sine refocussing pulse to cause a 

15 phase evolution of a plurality of nuclear spins in said patient; 

3. creating a second slice selection magnetic field gradient 
simultaneously with application of the 180° RF sine 
refocussing pulse to cause phase reversal only within a 
second slice of tissue within said patient intersecting the first 

20 slice, thereby creating a rod-shaped region of excited nuclear 

spins with reversed phase, 
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4. applying a second set of diffusion gradients in the 'X\ "Y' 
and "Z" directions after application of the 180° RF sine 
refocussing pulse to cause a second phase evolution of the 

25 nuclear spins, thereby encoding diffusion of the nuclear 

spins in the excitation region, and 

5. applying a readout gradient to said rod-shaped excitation 
region along the length of the excitation region being the "Z" 
direction for which a temperature-sensitive profile is desired; 

30 c) receiver means for receiving an MR response signal from the 
rod-shaped excitation region; and 
d) computation means for computing a temperature vs. position 
profile along the readout gradient from the MR response 
signal received from the receiver means. 

4. The MR pulsed heat system of claim 1 wherein the MR imaging 
means comprises: 
a) a radiofrequency (RF) transmitter for transmitting an RF 
excitation pulse and an RF refocussing pulse into said patient; 
5 b) a gradient means for. 

1. producing a time-varying magnetic field gradient in an "X" 
direction and time varying magnetic field gradient in a "Y" 
direction orthogonal to the 'X* direction, both gradients 
applied simultaneously with the RF excitation pulse exciting 

10 longitudinal magnetization of an elongated excitation region 

in said patient, 

2. applying a first set of diffusion diffusion gradients in the "X", 
T directions and a "Z w direction orthogonal to both the "X M 
and "Y" directions after excitation of the excitation region 

15 and before application of the refocussing pulse, 

3. applying a second set of diffusion gradients in the X, "Y" 
and "Z" directions after application of the first set of diffusion 
gradients, and 

4. applying a readout gradient to the elongated excitation region 
20 along a direction for which a temperature-sensitive profile is 

desired; 
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c) receiver means for receiving an MR response signal from the 
elongated excitation region; and 

d) computation means for computing a temperature vs. position 
25 profile along the readout gradient from the MR response 

signal from the receiver means. 

5. The MR pulsed heat system system of claim 1 wherein the 
pulsed heat-producing means comprises: 

a) a source of pulsed optical energy; 

b) an invasive device adapted for insertion into the patient to 
5 reach the specific tissue of the patient; and 

c) an optical fiber having an external end and an internal end 
adapted to be fitted into the invasive device with its internal end 
within the specific tissue, the fiber being adapted to pass the 
optical energy from the pulsed optical energy source into the 

10 external end, to its application point in the specific tissue 
adjacent the internal end. 

6. The MR pulsed heat system system of claim 5 wherein the 
source of pulsed optical energy comprises a laser. 

7. The MR pulsed heat system system of claim 1 wherein the 
pulsed heat-producing means comprises: 

an ultrasonic transducer adapted for generating pulsed 
ultrasonic energy concentrated at a focal point, the focal point 
5 being the application point. 

8. A method of performing heat surgery on a patient, as guided by 
magnetic resonance (MR) imaging, comprising the steps of: 

a) creating an internal image of tissues of said patient; 

b) determining the position of a specified tissue in the patient; 

5 c) applying pulsed heat at a predetermined pulse frequency to the 
specified tissue to create a heated region within the specified 
tissue; 
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d) monitoring the heated region with a magnetic resonance (MR) 
imaging means to acquire real-time temperature-sensitive 

10 images; and 

e) adjusting the frequency and the position of the heated region so 
as to heat the specified tissue without substantial injury to 
adjacent tissue. 

9. The method of performing heat surgery of claim 8 wherein the 
step of monitoring the heated region comprises the steps of: 

a) transmitting an RF excitation pulse into said patient; 

b) producing a first time-varying magnetic field gradient in an 
5 "X" direction and a second time-varying magnetic field 

gradient in a "Y" direction orthogonal to the "X" direction, 
both gradients applied simultaneously with the KF excitation 
pulse so as to excite transverse magnetization in an elongated 
excitation region in said patient; 
10 c) applying a dephasing magnetic field gradient to said patient 
for a predetermined duration having a predetermined 
amplitude of a first polarity, and 

d) applying a readout magnetic field gradient in a "Z" direction 
along the length of the elongated excitation region, having a 

15 second polarity opposite the first polarity, and having a product 
of duration multiplied by amplitude substantially larger than a 
product of duration multiplied by amplitude of the dephasing 
magnetic field gradient; 

e) receiving simultaneously with the application of the readout 
20 gradient, an MR response signal M^g which is equal to 

longitudinal magnetization Mz, defined by: 

M 0 for the first MR response signal 

for subsequent MR response signals 



M z = 



Mo 



l-e T > 



from the elongated excitation region, where Mo is the 
equilibrium magnetization, Ti is the spin lattice relaxation 
25 time, and Tr is the pulse-sequence repetition time; 
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f) waiting a time period Tr, as measured from the beginning of 
the transmission of the RF excitation pulse; 

g) repeating steps "a" through T for at least one repetition; 

h) computing a temperature-sensitive profile along the readout 
30 gradient, employing the MR response signal M a i g , for any but 

the first MR response signal . 

10. The method of performing heat surgery of claim 8 wherein the 
step of monitoring the heated region comprises the steps of: 
a) transmitting an RF sine(x)/x "sine" selection pulse into said 
patient; 

5 b) applying a slice selection magnetic field gradient along an "X" 
direction sinrultaneously with application of the RF sine 
selection pulse to cause excitation of a first slice of tissue 
within said patient; 

c) applying a first set of diffusion gradients in the "X" direction, 
10 and a T and "Z" direction orthogonal to the T direction to 

cause a phase evolution of a plurality of nuclear spins in the 
first slice of said patient; 

d) transmitting a 180° RF sine refocussing pulse into said patient 
reversing the phase of the nuclear spins; 

15 e) applying a second slice selection magnetic field gradient 
simultaneously with the application of the 180° RF sine 
refocussing pulse causing excitation of a second slice of tissue 
within said patient intersecting the first slice and thereby 
create a rod-shaped excitation region; 

20 f) applying a second set of diffusion gradients in the T, T and 
"Z" directions after application of the 180° RF sine refocussing 
pulse to cause a second phase evolution of the nuclear spins 
and thereby encode diffusion of the nuclear spins; 
g) applying a readout gradient to the rod-shaped excitation region 

25 along the length of the excitation region, the rod-shaped 
excitation region being along the "Z" direction for which a 
temperature-sensitive profile is desired; 
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h) receiving an MR response signal from the rod-shaped 
excitation region; and 
30 i) computing a temperature vs. position profile along the readout 
gradient from the MR response signal. 

11. The method of performing heat surgery of claim 8 wherein the 
step of monitoring the heated region comprises the steps of: 

a) transmitting an RP excitation pulse into said patient; 

b) applying a time-varying magnetic field gradient in an T 

5 direction and a time-varying magnetic field gradient inaT 
direction orthogonal to the "X" direction, both gradients 
applied simultaneously with the RF excitation pulse so as to 
excite longitudinal magnetization of an elongated excitation 
region in said patient; 
10 c) applying to said patient a first set of diffusion gradients in the 
T and HT directions and a M Z W direction orthogonal to both 
the T and "Y" directions after excitation of the excitation 
region; 

d) transmitting an RF refocussing pulse into said patient; 

15 e) applying a second set of diffusion gradients in the "X", T and 
"Z" directions after application of the RF refocussing pulse; 

e) applying a readout gradient to said elongated excitation region 
in a direction which a temperature-sensitive profile is desired; 

f) receiving an MR response signal from the elongated excitation 
20 region; and 

g) computing a temperature vs. position profile along the readout 
gradient from the MR response signal. 

12. The method of performing heat surgery on a patient of claim 8 
wherein the step of applying pulsed heat comprises the steps of: 

a) focus sing ultrasonic waves on an application point; and 

b) positioning the application point at a predetermined location 
S within the specific tissue. 
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13. The method of performing heat surgery on a patient of claim 8 
wherein the step of applying pulsed heat comprises the steps of: 

a) creating optical energy; and 

b) passing the optical energy through an optical fiber to a 

5 predetermined application point within the specific tissue. 
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